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Abstract 

We present simulation results of 7 + Jet analysis using CMS (Compact Muon 
Solenoid) Object-Oriented software at the Large Hadron Collider (LHC) center of 
mass energy ^/s=14 TeV. The study of direct photon production helps in validating 
the perturbative Quantum Chromodynamics (pQCD) and providing information 
on the gluon distribution in the nucleons. Direct photon processes also constitute 
a major background to several other Standard Model (SM) processes and signals 
of new physics. Thus these processes need to be understood precisely in the new 
energy regime. In this work, we have done a detailed study of the GEANT4 simu- 
lated 7 -t- jet events generated with Pythia, and the related background processes. 
Isolation cuts have been optimized for direct photon which improves the signal over 
background ratio by ~ 25% as compared to previous studies done in CMS. The 
inclusion of a large A0 cut between the photon and the leading jet at 40*^ in the 
analysis leads to a further increase of ~ 15% in S/B, thus giving an overall gain of 
~ 42% in S/B ratio. 
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Introduction 



Direct or prompt photons are photons that emerge directly from the hard scattering of the 
partons and not from secondary decays or as the radiation product of initial or final state 
partons. Since these photons come directly from parton-parton hard scattering, they provide 
a clean signature of the underlying hard scattering dynamics. Prompt photon production in 
hadronic interactions provides a precision test of pQCD predictions as well as information on 
the gluon density inside the colliding hadrons. Theoretical investigations have shown [1] that 
the cross-section for direct photon production has a pseudorapidity (r/) dependence which is 
sensitive to the parameterization of the gluon distributions[2j. The LHC energy will provide an 
opportunity to determine gluon density in a proton in new kinematic region of x{2 x 10~^ < 
X < 1.0) and (1.6 x 10^ - 2 x 10^{GeV/cf)^. The analysis of the LHC data in combination 
with results from the Tevatron and HERA would allow to extend the QCD analysis in region 
102 < g2 < io^(GeV/cf[A\. 

The study of physics beyond the SM at the LHC also requires a complete understanding 
of the QCD processes. Some of the QCD processes constitute major background to other SM 
processes, and also to signals of new physics. Thus these processes need to be well understood 
precisely in the new energy regime. In this work, we have concentrated on the QCD process 
pp — y 7'^*'^ + 1 jet + X, where X can be anything. This process is a major background to 
other Standard Model processes such as H — > 77 [S], W"^ — >• 7r^7[6], and signatures of physics 
beyond the SM such as large extra dimensions O [8] and SUSY[9]. 

Inclusive direct photon production has been studied extensively by various experiments 
from fixed target energies PH] to Tevatron collider [TTl [T21 [131 [H [IS]. A recent study by the DO 
experiment at Tevatron has reported on the measurment of the triple differential cross-section 
for pp — > 7'^*'' + jet at -^5=1. 96 TeV|16j. In this work, the photon + jet cross-sections has been 
compared to the next-to leading order (NLO) based on the program JETPH0X|17^ il8j with 
CTEQ6.1M PDF set of parton distribution functions|19[ [2U]. However, both the completeness 
of theoretical calculations and consistency of the available data have been subject of intense 
discussion. The results from various experiments have shown a definite pattern of deviation 
from theory at low Pt where the measured cross sections have a steeper slope. The origin of the 
disagreement has been attributed to the effect of initial-state soft gluon radiation which is not 
properly accounted for in the theoretical calculations. It has been shown that the inclusion of 
some additional transverse momentum "/c^" smearing of the partonic system due to soft-gluon 
emission yields better description of the data[2Tl [22]. These effects are however expected to be 
negligible in the kinematic range explored in the present analysis. 

At the leading order, the direct photon production is defined by two QCD processes; Quark- 
Gluon Compton scattering, q + g — > 7'^*'' + q and Quark-Antiquark annihilation, q + q — > 
^dtr _|_ g J^wQ to abundance of low momentum fraction gluons in the proton at the LHC 
energy, Compton scattering becomes the dominant process contributing to the prompt photon 
production over most of the kinematical region. 

Unfortunately, the advantages of photon as a clean probe of parton distributions are offset 
by large QCD backgrounds which are ~ 10^ to 10^ times larger than that of direct photon 
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signal. The background contribution to direct photon is mainly caused by the events where 
high Pt photons are produced in the decays of neutral mesons such as 7r°, t], and tu"; and 
from events where the photons are radiated from the quark (i.e, bremsstrahlung photons in the 
next-to-leading order QCD subprocesses such as qg — > qg, qq — > qq and qq — > qq). Any 
analysis must separate true direct photons (those coming from the hard scattering) from those 
copiously produced in the decays of 7r°, 7], /c°, uj^ and bremsstrahlung photons emerging from 
high Et jets. Isolation cuts imposed on reconstructed photon candidates effectively suppress 
these backgrounds |23j. The purpose of this work is to study the 7 + jet events and its back- 
ground processes with full detector simulation and reconstruction, and optimize the signal over 
background (S/B) ratio at the LHC energy. 

In this work, we present the GEANT4 simulated results of Level-3 Triggered photons (also 
called High Level Trigger (HLT) photons) using CMS Software packages. The event generation 
for 7'^*'"+ jet signal and background processes has been done using CMKIN[24j. The passage 
of particles through the detector geometry, showering, energy loss in the calorimeters and 
reconstruction of the events are modeled in the CMS simulation and reconstruction packages 
namely, OSCAR[25] and ORCA[26j respectively. The simulated result has been compared with 
the theoretical calculations. 

The rest of the paper is organised as follows. In the next section, we discuss the physics of 
7+ jet. Section 3 gives a brief description of the CMS detector relevant for the present analysis. 
In section 4, we discuss event generation. Section 5 describes photon isolation and its effect on 
signal and background. In section 6, we discuss the effect of A0 cut between the photon and 
the jet on S/B. Statistical and systematic uncertainties are discussed in section 7. In the last 
section, we summarize this analysis with the conclusions. 

The Physics of 7 + jet 

Fig.[T]shows the two leading order (LO) processes namely, "Quark- Antiquark Annihilation" and 
"Quark-Gluon Compton Scattering" which contributes to 7 + jet events. In quark-antiquark 
annihilation, a direct photon is produced along with an outgoing gluon while in quark-gluon 
Compton scattering, a gluon scatters from a quark producing a photon in the final state along 
with the scattered quark. An example of a next-to-leading order (NLO) real correction and 
corresponding virtual corrections (which cancel the soft and coUinear divergences arising from 
the real radiation) are shown in fig. [2J 

If the initial parton energies are equal, the photon and the jet are produced back to back 
(?7 ^ 0). In case of large imbalance between initial parton energies, the photon and the jet tend 
to be produced at small angles, either both forward or backward. Since the gluon momentum 
distribution increases monotonically at lower momentum fraction values [27], high r] photons are 
likely to come from gluon Compton scattering process rather than annihilation process. 

The main background to the direct photon signal comes from jet fluctuations. While most 
of the jets consist of many particles, thus easily distinguishable from a single photon, a small 
fraction (~ 10^'^ - 10"^) of jets fragments in such a way that a single particle carries most 
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Figure 1: Leading order processes contributing to direct photon production. 




Figure 2: NLO real diagram and corresponding virtual diagrams contributing to direct photon 
production. 



of the momentum of the parent parton. The lightest and therefore most commonly produced 
neutral meson, the 7r° decays into two photons with a branching ratio of ~ 99%. Photons 
produced in decays of high 7r° (and other neutral mesons) are very close to each other 
in the electromagnetic calorimeter forming a single shower thus mimicking a direct photon. 
Therefore, non-direct photons are common in hadronic collisions. Since jet production rates 
are ~ 10^ times larger than that of direct photons (depends on I77I and Py), the number of jets 
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which fragment to a single meson and then decay to photons, referred as electromagnetic (EM) 
jets, are comparable to the level of direct photons. The background to direct photon + jet can 
also come from photons produced in initial and final state radiation. In these cases, the photon 
is not produced directly from the interaction vertex, and therefore, is not really a direct photon. 
In order to reduce the background contribution, an isolation criterion is imposed on the photon: 
typically a cone of size R (in pseudorapidity {r]) and azimuthal angle (0)) around the photon 
is required to have less than a certain amount of energy. Since a photon produced through 
bremsstrahlung tends to be coUinear with the parent quark (and thus the jet resulting from 
the quark's subsequent fragmentation), the requirement that the photon be isolated removes 
all but very large angle radiations. Fig. [3] and H] show major backgrounds to direct photons, 
namely, the EM jets and the bremsstrahlung. 



q 




Figure 3: Background contributions to Direct Photon Events from EM jets. 



The CMS Detector 

The Compact Muon Solenoid (CMS) is a general purpose detector with in steridian cover- 
age for efficient detection and precise energy /momentum measurement of electrons, photons, 
muons, hadrons (jets) and missing transverse energy. From inside to outside, it comprises of 
a pixel detector, a silicon microstrip tracker system, a lead tungstate crystal electromagnetic 
calorimeter (ECAL), a sampling hadronic calorimeter (HCAL), a 4 Tesla soleniodal magnet 
and the muon detectors. 

The innermost component of the tracking system is the highly segmented silicon pixel detec- 
tor which consists of 3 layers of barrel complemented with 2 disks on each side in the endcaps. 
It consists of 66 million pixels of size 100 /xm x 150 /xm yielding a three dimensional spatial 
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Figure 4: Background contributions to Direct Photon Events from final-state radiation pro- 
cesses. 



resolution of ~10 /im in the r-0 and ~20 /im in the z-coordinate, which will play a crucial 
role in the vertex identification and track reconstruction. The following part of the tracking 
system consists of silicon microstrip detector modules arranged in ten concentric layers in barrel 
and nine disks in the endcaps providing a single point resolution between 23-52 /xm in r-0 and 
230-530 /im in z. It provides the necessary lever arm for accurate momentum measurement 
and improves the precision of vertex measurement. The tracker has geometrical coverage upto 
\ri\ = 2.5. The momentum resolution APt/Pt of high Pt charged particles in the tracker is 
~ {15{PT/TeV) + 0.5)% for \ri\ < 1.6 and becomes ~ (60(Px'/TeV) + 0.5)% as |?7| approaches 
2.5[28l[^|30]. 

The electromagnetic calorimeter is a hermetic, homogeneous calorimeter consisting of 61,200 
lead tungstate {PbWOi) crystals mounted in central part corresponding to < \ri\ < 1.479, 
and 7324 crystals mounted in each of the two endcaps in the range of 1.479 < \ri\ < 3.0. 
The crystals are 23 cm long, have short radiation length(Xo =0.89 cm) of Moliere radius(21.9 
mm) and correspond to ~26Xo, thereby containing more than 99% of the shower energy. The 
granularity of the calorimeter is At] x A0 = 0.0175 x 0.0175 in the barrel and upto Arj x 
A(f) = 0.05 X 0.05 in the endcaps. The resolution for single photons with energies of 25- 
500 GeV can be parameterized as as/E = 3.6%/\/E ® 124/ E © 0.26%. An additional Pb-Si 
preshower detector would be installed in front of the endcaps covering 1.65 < |?7| < 2.61 for 
efficient rejection of -7 /tt" [301I3T] . 

The sampling hadronic calorimeter is made up of plastic scintillator tiles inserted between 
brass absorber plates and extends upto |?7| < 3.0. The barrel (endcaps) consists of 5 cm (8 cm) 
thick brass plates interleaved with 3.8 mm thick scintillator tiles. The barrel calorimeter has 
a depth of 79 cm corresponding to 5.15 hadronic interaction lengths (A/) and a segmentation 



5 



of At] X A0 = 0.087 X 0.087. Outside of the soleniod, an additional layer of 10 mm thick 
scintillator called outer hadron calorimeter covering the \ri\ < 1.26 region, is placed to ensure 
sampling of hadronic shower with ~ 11 A/. Hermecity of the HCAL is extended to |?7| < 5.0 by 
adding forward detector made up of steel absorbers and embedded radiation hard quartz fibers, 
which provide a fast collection of Cerenkov light. The resolution of the HCAL ranges from 
cte/E = 65%/a/E ©5% in the barrel (at r] = 0) and cte/E = 83%/VE © 5% in the endcaps, 
to as/E = 100% /y/E © 5% in the forward region [301 E2]. 

The muon detector consists of drift tubes in the barrel {\ri\ < 1.2) and cathode strip chambers 
in the endcaps (1.2 < \ri\ < 2.4) for precision track measurement. A set of resistive plate 
chambers are used in both parts of the system for resolving individual bunch crossing time. The 
spatial resolution of the muon system varies between 50 to 200 fim and standalone momentum 
resolution is at most 15% for a 10 GeV Pt muon and 40% at 1 TeV[30l 133]. 

Event Generation 

At the LHC center of mass energy, the cross-sections for the channels pp — > 7+ jet and pp 
— > jets are very large, thus making it nearly impossible to simulate direct photon signal and 
its background processes equivalent to an integrated luminosity of ~ lfb~^. In order to achieve 
this, a selection based on kinematic information of the particle is performed at the generator 
level which allows only those events to be fully simulated that are likely to pass the analysis 
criteria. In our study we have used ItiemcalI < 2.6 to match the preshower coverage. The 
CMS experiment has optimized the LI trigger cut for a LI isolated single photon to be 23 
GeV, while that for the HLT photon it is 80 GeV at a luminosity of 2 x 10^^cm"^s~^[3n]. For 
the signal, a pre-selection has been used at the generator level which retains only those events 
having a direct photon with P^ > 70 GeV and \r]'^\ < 2.8, which are well below the final analysis 
cut of P^ > 80 GeV and \ri'^\ < 2.6. This preselection helps in discarding those events before 
simulation which are unlikely to pass geometrical acceptance cut. The LI trigger efficiency for 
the signal has been estimated to be above 99% over the entire kinematic range [30]. 

For the backgrounds, we have used Pythia generated QCD dijet events. The samples 
{eg03-jets-le) have been preselected at the generator level with an enhanced LI trigger 
efficiency [23]. In this pre-selection, the first step is to look for the seed particles of electro- 
magnetic objects like photons, electrons and positrons, which have Pt > 5 GeV and |?7| < 2.7. 
Candidate electromagnetic calorimeter trigger tower energies are then estimated by adding en- 
ergies of all electromagnetic particles found within a cone of At] <0.09 and A0 <0.09 from the 
seed. Trigger tower candidates that lie within Arj < 0.2 and A0 < 0.2 from each other are sup- 
pressed and only those with the highest Pt are retained. The LI single photon electromagnetic 
trigger is simulated by requiring that one such candidate has transverse energy larger than 20 
GeV. 

Signal and background events have been reconstructed using ORCA version 8.7.3. In the 
electromagnetic calorimeter, the Hybrid and Island clustering algorithms |34j were used to con- 
struct basic clusters in the barrel and the endcap regions respectively. A supercluster is a cluster 
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of the basic clusters. A photon candidate is associated with each reconstructed supercluster. 
Events were generated in two different Pt bins: 50 GeV < Pt < 170 GeV, and Pt > 170 GeV. 

The cross-sections, number of events and corresponding integrated luminosity for signal and 
background samples are shown in Table [H The samples have been generated using Pythia[35]. 
simulated and reconstructed using OSCAR and ORCA. We require that the leading photon 
should have Pj^ > 80 GeV and \ri'^\ < 2.6. Fig. [5] shows the number of events/GeV after P^ 
and 1] requirements calculated as a function of P-f! for the 7 + jet signal and its background for 
an integrated luminosity of 1 fb"^. 



Data 
Samples 


Pt 
(in GeV) 


a generated 
(m pb) 


a preselected 
(m pb) 


Number of 
Simulated 
Events 


Integrated 
Luminosity 

(in pb^^) 


7 + jet 


50-120 


8.53E+03 


1.82E+03 


30000 


16.5 


7 + jet 


120-170 


2.77E+02 


2.35E+02 


12740 


54.2 


7 + jet 


170-230 


6.76E+01 


6.19E+01 


11918 


192.5 


7 + jet 


230-300 


1.90E+01 


1.83E+01 


11531 


630.9 


7 + jet 


> 300 


9.07E+00 


8.98E+00 


9887 


1212.3 


eg03_jets_le_50170 


50-170 


2.43E+07 


4.35E+06 


2314853 


0.53 


eg03_jets_le_170up 


> 170 


1.34E+05 


1.07E+05 


473668 


4.4 



Table 1: Cross-sections, number of simulated events, and equivalent integrated luminosity for 
7 + jet events and its background samples. 



Photon Isolation 

As shown in fig. [5l jets misidentified as photons give a very large background to the direct 
photon signal. Isolation cuts have been optimized and applied to reject photons coming from 
jets by detecting other fragments of the jet which accompany the photon. 

Tracker Isolation 

Direct photons produced in hard scattering are not accompanied by any charged particles, 
whereas in cases where a jet fragments into a 7r° mimicing a photon or the jet bremms a 
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Figure 5: Number of Events/GeV for the 7 + jet signal and its backgrounds for an integrated 
luminosity of 1 fb~^. No isolation requirements have been applied on the photon. 



photon, the electromagnetic object is usually accompanied by tracks of low momentum charged 
hadrons. Thus the tracker isolation is an important parameter in differentiating the signal 
from the background events. The tracker isolation criteria is based on the number of charged 
particle tracks above a certain Pt threshold calculated in a cone R (R=v^Arp~+"A0^) around 
the photon candidate. The algorithm contains following three parameters: 

1. Cone size R: The size of the cone R around the photon candidate, in which the number 
of charged tracks is counted. 

2. Pt threshold (P-^^'"^'*): Only tracks greater than threshold Pt, P^^*"^*, are considered for 
photon isolation. 

3. Number of tracks threshold {N^^'^'^): If the number of tracks in cone R with track 
Pt > Pt^^^'^ is larger than Ntk"^^, then the photon candidate is considered non-isolated, 
otherwise it is considered as an isolated photon. 

We have done a detailed study of the effect of Pt threshold of the track (P^'*'"'^'*) and 
the number of allowed tracks [N^^'^^) inside the isolation cone (R), on signal efficiency and 
background rates. The study has been done separately for the barrel and the endcaps. The 
results for the barrel and the endcaps are shown in fig.[6]and fig.[7]respectively. Now limiting the 
discussion to fig. [6]^a) only, the six curves show the signal efficiency as a function of background 
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rate for three different values of track Pt threshold (1.0 GeV, 1.5 GeV and 2.0 GeV) for each 
of the two cone sizes (R=0.3 and R=0.4). From left to right, the five points on each of these 
curves show the change in signal efficiency and the background rate as the number of tracks 
allowed in the vicinity of the photon increases from to 4. The background rate is found to be 
very sensitive to the number of tracks allowed inside the cone. For example, if the number of 
allowed track changes from to 1, the signal efficiency increases by 4-10% while the background 
rate goes up by a factor of ~ 2 in the barrel and ~ 1.5 in the endcaps. Hence, for the maximum 
background rejection, N^^'^'^ is fixed to zero. For a fixed cone size, both signal efficiency and 
background rate increase with increase in track Pt threshold. Due to larger tracker inefficiencies 
the performance of the tracker variables is not as efficient in endcaps. Different values of P^^^'^'^ 
and cone size have been taken into account to estimate the optimized S/B ratio. 
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Figure 6: Performance of the tracker isolation variables for different track P^^'^'^^ in the barrel 
for cone sizes R =0.3 and 0.4 (a); R = 0.5 and 0.7 (b). 



Electromagnetic Calorimeter Isolation 

The ECAL isolation is based on the sum of transverse energies deposited in the basic clusters in 
cone of size R in — plane around a photon candidate. The basic clusters that belong to the 
photon candidate supercluster are not counted as a part of the sum. If the sum of transverse 
energies {T^Etecal) is below a certain threshold value {E^ECALly the photon candidate is 
considered isolated, otherwise it is considered non-isolated. 

The effect of making the values of -Etj^cal increasingly stringent on the signal efficiency 
as well as on background rate for different cone sizes R for the barrel and the endcaps can be 
seen in fig. |8] (a) and (b) respectively. From right to left, the value of E'^^q^j^ is gradually 
decreased from a very large value in steps of 0.5 GeV to yield the curves. Just to illustrate 
the nature of dependence of the signal efficiency and background rates, the plots shown have 
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Figure 7: Performance of the tracker isolation variables for different track P^^^"^^ in the endcaps 
for cone sizes R =0.3 and 0.4 (a); R = 0.5 and 0.7 (b). 
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Figure 8: Performance of the variable Etecal for various cone sizes R 
for the barrel (a); and for endcaps (b). 



=0.3, 0.4, 0.5 and 0.7 



variation in -Ej'^cal from 17 GeV to 1 GeV. The stringent selection cut on Y^Etecal provides 
a good background rejection. A tighter Etecal cut in the endcaps as compared to the barrel 
region can further enhance background rejection. 

Hadron Calorimeter Isolation 

Usually the hadronic activity around the background photons tends to be much more than 
around the direct photons. Thus, the HCAL isolation can also be used for photon selection. 
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For isolation in HCAL, the sum of transverse energies of all the particles depositing energy in 
a cone of size R around a photon candidate is calculated. If the sum of transverse energies 
of the hadronic particles, T^Ethcal is below a certain threshold value {E^h'^al)^ ^^e photon 
candidate is considered isolated. Since the information from the charged mesons and hadrons 
is used for tracker isolation, HCAL isolation offers partially redundant information. 

Fig. [9] (a) and (b) show the effect of variation in the E'^^^\j^ on the signal efficiency and 
the background rate for different cone sizes for the barrel and the endcaps respectively. From 
right to left the E'^^^\i^ is gradually decreased from a very large value in steps of 0.5 GeV 
to yield the curves. The HCAL isolation variable offers a better background rejection in the 
endcaps. 

Combination of the Detectors 

After studying the isolation parameters separately, the task is to combine their effects such that 
the S/B ratio is maximized while retaining a very high signal efficiency. This study has been 
done for the various detector combinations (for example: Tracker + HCAL, Tracker + ECAL, 
and Tracker + ECAL + HCAL) for different threshold parameters. We observe that a better 
reduction in the background rate without compromising the signal efficiency can be achieved 
by combining the isolation requirements from various detectors. Some combinations which are 
found to yield better signal to background ratio are shown in Table [H Table [3] shows the signal 
efficiency, signal and background rates, and S/B values for respective selection cuts mentioned 
in Table [21 It should be noted that the selection cut A has been used in the results mentioned 
in CMS Physics TDR-I[30]. With selection cut C, the background rate is further reduced by 
~ 22% with only 1% loss in the signal efficiency leading to a gain of ~ 26% in the S/B ratio. 
Fig. [TU] shows the number of events/GeV for an integrated luminosity of 1 fh~^ calculated as 



11 



a function of photon Pj- for the 7 + jet signal and its background for selection criteria C. The 
effect of the isolation requirement on the number of events/GeV can be observed by comparing 
fig. [5] and fig. [M 



Selection 
Cut 


Cone Size R 


TrackPf'^^'^^ 
(GeV) 


Tp'l'hres 
^TECAL 

(GeV) 


rp'l'hres 
^THCAL 

(in Barrel) (GeV) 


rpThres 
^THCAL 

(in Endcaps)(GeV) 


A 


0.3 


1.5 


1.5 


6.0 


4.0 


B 


0.3 


1.0 


1.5 


6.0 


5.0 


C 


0.4 


1.5 


2.0 


7.0 


5.0 


D 


0.4 


1.5 


2.0 


6.5 


5.0 


E 


0.5 


1.5 


2.5 


8.0 


6.0 



Table 2: Selection Cuts corresponding to the different values of isolation parameters threshold. 



Selection 


Signal 


Signal 


Background 


S/B 


Cut 


Efficiency 


Rate (Hz) 


Rate (Hz) 


Ratio 


A 


0.76 


2.12 


1.40 


1.52 


B 


0.76 


2.10 


1.26 


1.66 


C 


0.74 


2.09 


1.09 


1.92 


D 


0.73 


2.07 


1.06 


1.94 


E 


0.70 


1.99 


0.84 


2.37 



Table 3: Signal efficiency, signal and background rates for various selection cuts mentioned in 
Table [2] for low luminosity (L = 2 x 10^^ cm-^s"^) at ^^5=14 TeV. 

Fig. [TT] shows the the predictions for the differential cross-section of the 7 + jet events as a 
function of Pj^ compared with theoretical LO and NLO calculations [36] • The error bars on the x- 
axis represents the bin width. The differential cross-section calculated from the simulated data is 
represented as LO (calculated after full simulation) in the above mentioned figure. To evaluate 
it, the cross-section from Pythia was generated with CTEQ5L parton distribution function 
and renormalization scale n = Pt- Various factors like geometrical acceptance, efficiency of 
the analysis cuts and luminosity have been taken into account. Theoretical LO and NLO 
calculations are provided by J. Owens [37]. The results are found to be in good agreement with 
the theoretical LO calculation. As expected the NLO QCD contribution is higher than the LO 
in the whole Pt range under analysis. 

Acj) Cut 

In 7 -|- jets events, the jets have been reconstructed using Iterative Cone type algorithm with 
cone size R =0.5[38]. Since the HCAL extends up to |?7| < 5.0 the selection cuts used on the 
highest Pt jet are : P^"^ > 40 GeV and {rj^^^l < 5.0. Usually 7 -|- 1 jet events are produced in a 
"back-to-back" fashion, but this topology is disturbed by initial-state and final-state radiation. 
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The events with the vector Pj,^ being "back-to-back" to the vector Pj. within A0 can be defined 
by the equation; 

$(7, jet) = 180° ± A</) 

where $(7, jet) is the angle between Pj^ and Pj^^* vectors. To further optimize S/B ratio while 
retaining a very large signal efficiency, the effect of A0 cut has been studied after applying 
all photon isolation requirements. Fig. [T^ shows the signal efficiency as a function of A0 cut. 
Fig. [13] shows the performance of A0 as a function of signal efficiency vs. S/B ratio for all 
selection cuts noted in Table [2l when A0 is varied in steps of 5°. It is found that the limit 
A(f) < 40° can be used as an additional cut to improve the S/B ratio. Requiring A0 < 40° 
improves the S/B ratio by ~ 15% with a bare reduction of 1 — 2% in signal efficiency. 




|-i . I I .... I 1 , !■ , ■! I: , , , , 1 , ,11,1 Nil 1,1,11 1 III I II I, II, 1 , , ,J. 

100 200 300 400 500 600 700 800 900 1000 



Figure 10: Number of Events/GeV for the 7 + jet signal and its background for an integrated 
luminosity 1 fb^^ after applying the selection cut C. 



Statistical and Systematic uncertainties 

The statistical uncertainty on the signal and the background rates are found to be ~ 1% 
and ~ 5% respectively. We have studied some of the potential sources of theoretical and 
experimental uncertainty. The uncertainty on the signal (background) rate arising from the 
choice of the pdf used in the event generation is found to be ~ 3% (~ 5%). The variation of 
the hard scale from default choice of in Pythia to Q"^ = s (MSTP 32=4) results in ~ 7% 
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Figure 11: The LO (after full simulation) and theoretical LO and NLO calculations for the cross- 
section of 7 + jet events as a function of Pj-. The calculated points have been complemented 
by lines fitting the corresponding cross sections. 



uncertainty on the signal rate and ~ 27% on the background rate. The estimated uncertainty 
due to photon trigger is ~ 1%[3 EH]- For Pj^^ > 60 GeV, the uncertainty due to jet energy 
scale is ~ ±7%. Experimental measurement uncertainties such as photon energy scale and jet 
energy resolution effect the rates within ~ 1%. The systematic uncertainty for the signal and 
the background rate arising from the preselection efficiency is less than 3%. These uncertainties 
when added in quadrature leads to ~ 11% uncertainty on the signal rate and ~ 29% on the 
background rate. 

In the initial phase of the LHC operation, for an integrated luminosity of lfb~^, the error 
on the measured luminosity at CMS is expected to be ~ 10%[in]. Studies from other colhder 
experiments im IT5| ST] have reported uncertainties due to photon purity, inefficiency in event 
vertex determination and conversion of photons in the detector. However, the understanding 
of such sources of systematic uncertainty in the context of this analysis would be better and 
clear after the availablity of the LHC data. 

Conclusions 

In this analysis, we have explored the direct photon + jet production at the CMS in the region 
of photon transverse momentum (Pj^) from HLT Pt threshold of 80 GeV to a few hundred GeV. 
We have studied the tracker, the ECAL and the HCAL based isolation conditions for the direct 
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Figure 12: Signal efficiency as a function of variable for various selection cuts mentioned in 
Table [2 



photon physics. For low luminosity phase of the LHC operation, it is found that with optimized 
isolation conditions the background can be reduced by two orders of magnitude while the signal 
efficiency remains between 70% — 80%. Isolation requirement within a larger cone around a 
photon provides a clean and pure signal. It is found that for the entire Pt range in the analysis, 
carefully chosen values of the isolation threshold parameters can yield appreciable increase in 
the S/B ratio as compared to previous studies while retaining very high signal efficiency. We 
have found selection cut C (Table [2]) to be optimum. We have matched our results from the 
Pythia based simulation with an independent theoretical calculation and they are in excellent 
agreement. The topology of the 7 + jet events is such that at high Pt most of the direct photon 
and jet events are "back to back" within certain A0. The inclusion of a very wide A0 cut at 
40° in the analysis leads to a further increase of ~ 15% thus leading to an overall gain ~ 42% 
in S/B ratio. A better understanding of fake photon rejection will naturally lead to improved 
efficiency for detecting low mass Higgs decaying to photon pairs. 

Acknowledgments 

We are grateful to J. Owens for providing theoretical calculations. We thank M. Fieri for 
helping with the background samples. We also wish to express thanks to F. Demin, D. Futyan, 
A. Kumar, M. Mangano, C. Seez and R. Y. Zhu. The work has been supported by a grant from 
the Department of Science and Technology, Govt, of India. We gratefully acknowledge the 



15 



o 0.76 
0) 

"o 

0.74 

LU 

c 0.72 
0.7 



0.68 
0.66 



1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 

S/B 
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